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Abstract 

The generating functional of correlation functions for the XX Z spin chain is considered in 
the thermodynamic limit. We derive a system of integro-difference equations that prescribe 
this functional. On the basis of this system we establish the operator- valued Riemann-Hilbert 
problem for correlation functions of the XX Z spin chain. 

1 Introduction 

It has been shown by Korepin et al that many important problems arising from analysis of 
correlation functions of quantum solvable models are reduced to classical inverse scattering 
problems |lj. In particular, for the XXX spin chain, the generating functional of correlation 
functions is represented by the Fredholm determinants [0|, and a special correlation function, the 
so-called Ferromagnetic-String-Formation-Probability (FSFP), is proved to be connected with 
an operator- valued Riemann-Hilbert problem . With the help of the solutions of this problem 
one can compute the large distance asymptotic form of the FSFP Q. Physically, this makes 
clear the probability of finding a ferromagnetic string of adjacent parallel spins for a given value 
of the magnetic field in the XXX spin chain. 

In this paper, we consider the spin-g XXZ Heisenberg chain in the thermodynamic limit. 
The Hamiltonian is defined by 



HxXZ = ^On<+l + ^n+l + cos 2 Wn<+l 



hat 



:i.n 



where cr*i a n an d ^ are the Pauli matrices acting on the n-th site and h is an external magnetic 
field. The anisotropy cos 2n implies the critical regime. Any correlation function of the model 
can be obtained by means of the generating functional Q^ m '(a) ||. For example, a two-point 
correlation function {cr^af) is given by 



(VmVl) 
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(1.2) 
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with the lattice Laplacian defined by A m /(m) = f(m) — 2f(m — 1) + f(m — 2). This generating 
functional is represented by the Fredholm determinants: 

, ( ^|Mi-^»)y (L3) 

(vac | det(l - ^Vy m )\ a= o)\vac) 



Here the kernel is of the form: 



sin2r/ /exp(a - (^ 3 (A) + <^ 4 (^))e 1 (A)ei(^) + l 



sinh(A — n) \ sinh(A — /j, + 2ir/) 

_ e^" 1 (A)e 2 (^) + exp(q - y 3 (A) + y^G")) 
sinh(/x — A + 2i7/) 



;i-4) 



From now on the index m is omitted if unnecessary. The integration contour that appears in 
the Fredholm determinants lies on the real axis [—A, A], where a boundary value A is called 
the Fermi energy and depends on the anisotropy parameter 77 and the magnetic field h. The 
operators ei(A) and e2(A) are defined by 



m 



sinh(A + ir?)^ ^ i(A) / sinh(A + irj) \ (A) 



The operators ipj(X) (j = 1, . . . ,4) are bosonic quantum fields called the dual fields. They are 
decomposed into the momentum and coordinate fields: 

<p j (X)=p j (X)+q j (X) (j = l,...,4) (1.6) 

whose commutation relations are 

bi(A),Pfc(/x)] = fe-(A),gfc(/i)]=0 

bi(A),«*(M)] = U jk h(X,fx) + U kj h(n,X) (j,k = l,...,A) (1-7) 



sinh(A - /i + 2i7?) 

/i(A,ii) = log ^— . (li 

1 sin lr\ 



with 
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Notice that the dual fields are commutative: [fj(X),(pk(fJ>)] = (j, k = 1,...,4). They only 
produce vacuum expectation values according to the vacuum states defined by 

(vac I vac) = 1 (vac|<£j(A) = pj(X)\v&c) = 0. (j = l,...,4) (1-9) 

A special correlation function FSFP for the XXZ spin chain is defined by P{m) = 
Q( m '(— 00), where m gives the length of a ferromagnetic string [||, |6|. In this case the ex- 
pression of kernel (|l.4| ) is reduced to a simple form and enables us to formulate the related 
operator- valued Riemann-Hilbert problem in the same way as for the XXX spin chain |J. 
Furthermore, in the limit of strong magnetic field h — > h c = 4 cos 2 77, the kernel is related with 
the r-function of the Painleve V equation 0. On the basis of this fact we can evaluate any 



correlation function of the XXZ spin chain under an external strong magnetic field. Similar 
results hold in the asymmetric XXZ spin chain that is a non-hermitian generalization of the 
XXZ spin chain @, |. 

The aim of the paper is to derive a system of integro-difference equations that prescribe the 
generating functional Q( m > (a) and to establish the associated operator- valued Riemann-Hilbert 
problem. We do not take either the FSFP limit (a — ► — oo) or the strong magnetic field limit 
(h — > h c ), and treat directly the kernel ( |1.4| ) that contains the dual fields with four-species. This 
kernel is not of the desired type, however, it can be transformed into an integral operator that 



satisfies the integrable condition (see (2.7)). Thus the computation of the generating functional 
can be reduced to the solution of operator-valued Riemann-Hilbert problems. In section we 



rewrite the kernel (1.4) and introduce some integral operators. In particular the resolvent is 
expressed by products of two vectors (see Lemma ||) . In section || a system of integro-difference 
equations is derived. On the basis of this system we establish the operator-valued Riemann- 
Hilbert problem associated with the generating functional ( |1.3j ) in section |||. Section [|, the last 
section, is devoted to concluding remarks. 



2 Integral operators 



In this section we transform the kernel (1.4) into an integral operator that satisfies the integrable 
condition ( |2.7D and introduce some useful integral operators. 



Definition 1. We rewrite the kernel for the Fredholm determinant representation of the gen- 
erating functional (|1.3|) as 



V(x,y) 



d.s 
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aj(x\s)bj( 



(2.1) 



where vectors a(x\s) and b(x\s) are defined by 



a(x\s) 



\{q-q ) 



1 qx 



1/2 



bHx\s) 



i(q - q 



-uQ x 



1/2 



V 

/ 

V 



X 



mjl 



qx — 1 



i ■ qx — 1 

exp(-ig^s 



\ 



x -m/2 exp (_i q g2_± s +a + (pi ( x )_ (p 3 ( X )) 
x -m/2 exp ( [q -l<E_l s 



<f>2(x)) 



x—q 



:r "'"exp(iq ' ' ,. ■; i 



l qx-l 

x — q 
m/2 eM _ iq3 ^l s + ip2{x)) 

m/2 exp(-igf^ S + ^ 4 (x)) 



-m/2 

-x^^exp^- 1 ^-^ 
x 

X 



(2.2) 



(2.3) 



The superscript t indicates the transposition of a vector. The integration contour for variables 
x,y runs on the unit circle anti-clockwise: C : x = e l8 (—rp < 9 < 2tt + tp). The endpoints are 



denoted by (,± = e • The dual fields still obey the commutation relations (1/7) with 

-i 



h(x, y) = log ■ 



qx-l y-q 
x-q qy-1 



x-g gy-1 

qx-l y-q 



-1 



(2.4) 



The original kernel (|1.4|) is obtained by using transformations 

sinh(A + irj) 2 m ■ i sinh(A — irj) 



q = e *v V = i log — ttt rr (2.5) 



sinh(A — \rj) ' sinh(A + irj) 

and an identity that is valid for ^ < 77 < ir: 



usexp — 10 s + 10 s = i r. 2.5) 

Here an extra factor that remains the Fredholm determinants unchanged appear, however, it 
can be ignored because it has no influence on the representation of the generating functional 
(|1.3| ), We remark that the vectors a(x\s) and b(x\s) satisfy 

4 

dsJ2a j (x\s)b j (x\s) = 0. (2.7) 

i=i 

This relation is called the integrable condition and assures that the kernel ( |2.1| ) is well-defined for 
any x,y € C By virtue of this condition we can formulate the operator-valued Riemann-Hilbert 
problem in section |j. 

Hereafter some matrices that contain integral operators will appear. We call them the matrix 
integral operators and denote them by bold letters, except that "1" indicates the product of the 
delta-function and the unit matrix. The following convention for matrix integral operators is 
adopted unless it invites confusion: 

r /'OO 

(EF) jk (x, y\s, t) = dz / dr S^E jn (x,z\s,r)F nk {z,y\r,t). (2.8) 

JC Jo 

Definition 2. We introduce the vectors a(x\s) and (3(x\s) via 

((1 - ±V)a)(x\s) = a(x\s) ((3(1 - ±V))(x\s) = b(x\s) (2.9) 

and define the resolvent R(x,y) through 

(1 - iF)(l + &R) = (1 + &R)(1 - ±V) = 1. (2.10) 

Note that one of these relations can be derived by letting (1 + ^=R) act on both sides of the 
other relation from the left or the right. Thus (1 — -^V) and (1 + ^R) are commutative. 

Lemma 1. The resolvent is represented by 

r°° ds 4 

R(x,y)=i yWxIs^yls). (2.11) 

Jo x ~yfz{ 



Proof. The defining relation of the resolvent is rewritten as (1 — -^V)R = V. By using the 
identity (x — y) = (x — z) + (z — y) this relation is deformed into 

(x-y)R(x,y) - — dz dsJ~)aj(x\s)bj(z\s)R(z,y) 

27T Jc Jo ■—{ 

- -^ dzV(x,z)(z-y)R(z,y) 

= i / dsy>.,(2|s)6,(y|s). (2.12) 



Move the second term in the left side to the right side and recall the defining relations ( |2.9D 
( 2.10J ). The representation ( |2.11 ) is thus obtained. D 



Definition 3. We define the 4x4 matrix integral operators with variables s,t by 

l dx 
A jk (s,t) = i / —aj(x\s)b k (x\t) (2.13) 

Jc x 

f dx 

B jk (s,t) = i / — aj (x\s)f3 k (x\t). (j,k = l,...,4) (2.14) 

Jc x 

Lemma 2. The integral operator B(s,t) is the resolvent of A(s,t). Namely, 

(1 - ii)(l + iB) = (1 + iB)(l - iA) = 1. (2.15) 

Proof. The action of (1 — ^V) on -a(x\s) is computed as 



-a(x\s)-— / — V(x,y)a(y\s) 
x 2vr J c y 

= -a(x\s) + — / dy ( ) V(x,y)a(y\t 

x 2ir J c \x y) 



l((l-^A)a){x\s). (2.16) 



By virtue of this relation we obtain 

/ (IT* 

A jk (s,t) = i / — aj (x\s)(p k (l-±V))(x\t) 
Jc x 



c 

i / dx(5 k {x\t) ( ~aj(x\s) - — - / — V(x,y)a>j(y\s] 
Jc \ x 2n Jc V 



= i [ ^p k ( x \t)((l-±A)a) 3 (x\s) 
Jc x 

= «l-±A)B) jk (s,t). (j,k = l,...,4) (2.17) 

This means (1 — ^A)B = A and therefore implies the lemma. The proof has been completed. 

□ 



Applying Lemma |2| to ( p,16| ) we can obtain another representation of -a(x\s): 

la(x\s) = ~((1 + &B)a)(z\a) - ±- [ ^V(x,y)((l + ±B)a)(y\s). (2.18) 

X X Z7T J Q y 

We remark that the kernel V(x,y) commutes with A(s,t) and B(s,t) because variables are 
different. This relation plays an important role in the next section. 

3 Integro-difference equations 

In order to make clear the properties of the unknown vectors a( m \x\s) and (3^ m \x\s) we compute 
their dependences on m and ip, and derive a system of integro-difference equations for them. 
This helps us to formulate an operator-valued Riemann-Hilbert problem for the XXZ spin 
chain. 

Lemma 3. The vectors a.^ m '{x\s) and f3^ m '{x\s) obey the following integro-difference equations 
with respect to m: 

a ( m+1 \x\s) = ((^ 7l + ^(l-XA( m+1 )) 72 (l + ^B( m )))a( m ))(x| S ) (3.1) 

y/X 

(3^+ l ){x\s) = {p^)(^ l2 + l={l-^.A^) 11 {l + ^B^ + ^))){x\s) (3.2) 

where 71 and 72 are 4x4 diagonal matrices: 

7i = ^(l + ^ 2 ^) 12 = \{1-<J Z ®<J Z )- (3.3) 

Proof. By definitions of a^ m '(x\s) and b^ m '(x\s) their dependences on m are 

a ( m+1 \x\s) = (^ ll + ^- l2 ) a { - m \x\s) (3.4) 

y/X 

b( m+1 \x\s) = b( m \x\s)(Vx~ l2 + -^li). (3.5) 

y/x 

Similarly the kernel V^ m+1 '(x,y) is connected with V^ m '(x,y) as 

^ m+1 ) (x, V) = J^ (v (m) (x, y^^^ ^ b(m) (V\ s )^a {m) (x\s)\ . (3.6) 

Using these recursion relations we obtain 

a( rn+1 \x\s) - a( rn+1 \x\s) 
= _L(y( m+1 W m+1 ))(x|s) 

= # / %V^Xx,y)cc^Xy\s) - * (A^W" )^)- (3-7) 



It thus follows that 



x 2vr J c y/y 

= -L a ( m+1 \x\s) - — (A {m+1 ^ 2 a^)(x\s) 

yJX 2lTX 

= 7 ia (m) (*N) + ^((1 - ±A( m+ V) l2 a^){x\s) 

= ((1 - ^V^)W m) )(x\s) + |((1 " ^ F A(-+ 1 )) 72 (1 + ±BW) a W)(x\s) 

"^/ c f y(m)(:E ' 2/)((1 "^ A(m+1))72(1 + ^ jB(m))a(m))(y|s) - (3 ' 8) 

In the last equality the relation fl2.18|) is used. Thus the integro-difference equation (|3.1| ) is 
obtained by removing the action of (1 — ^V^ m '). In the same way ([T^) can be derived. D 

Lemma 4. The derivative of the vectors a(x\s) and f3(x\s) with respect to ip are 

d f a(x\s) = — J" \-^-r a (Ce\a) r dtJ2<xMt)Pj((e\t)) (3.9) 

d^/3(x\s) = -^J^l^^PiMpdt^ajiCelWMt))- ( 3 - 10 ) 



Proof. We show the differential equation (3.9). The other equation ( |3.1C| ) can be derived in the 
same way. According to the identity 



i^ f dxf(x) = J2 Cef(Cs) (3.H) 

J C L 



the differentiation of fl2.9| ) is computed as 

d^a(x\s) ~^~Y1 & V ( X > CeMCe\s) - ^{Vd^cx){x\s) = 0. (3.12) 



£=± 



By using the resolvent this is reduced to 

d$a(x\s) = — Y^ CeR(x, (e)ac(( E \s). (3.13) 



2vr 

£=± 



Because of Lemma [I] we obtain (|3.9|). □ 

The following lemma is not directly related with the formulation of operator-valued 
Riemann-Hilbert problems but useful for evaluating the large m asymptotic form of the gener- 
ating functional Q^ n '{a): 



Lemma 5. The Fredholm determinant det(l — ^V^™') satisfies the following recursion relation 
with respect to m: 

det(i-^y( m+1 ); 



exptrlog(l + ^ 72J B( m ); 



det(l-^yM) 
where the trace for matrix integral operators is defined by 



(3.14) 



/OO n 
]]ds k Yl K jun (s 1 ,s 2 )---K jnJl (s n ,s 1 ). (3.15) 



Proof. Due to the relation (|3.q) it follows that 

= 5{x-y)- — fvW(x,y)-ij -((3^(1 - ±V^)){y\s) l2 a^\x\s\ 
= ((l + ±G^)(l-±V^))(x,y) 
where 

r°° c\ s 

G^ m \x,y) = i / -f3^\y\s) l2 a^\x\s). 
Jo x 

Based on this recursion relation the Fredholm determinants of both sides are given by 

det(l - ^ (m+1) ) = det(l + ±G^) det(l - ^V™). (3.18) 



(3.16) 



(3.17) 



Here we note that \Jyjx in the left side has no influence on the Fredholm determinant. Take 
the logarithm of det(l + ^G): 



logdet(l + ±G) = trlog(l + &G) = £ ^yT^G"). 



(3.19) 



The trace of G n can be expressed in terms of B: 

tr(G ., = i-/nf f ft^ 

x 5Z a ji( 2; il' s i)/ 3 ji( a; 2|si)aj 2 (^2|s2)/5j 2 (x3|s2)---a : , n (x„|s n )/3 Jn (xi|s n ) 

jl,—,jn=2,3 



]^[ds fc ^ 5 ii,jn('Sl,Sn)5j 2 j 1 (s2,Sl)----Bi„,i n „ 1 (Sn,Sn- 



fe=l ji,...,j„=2,3 



= tr(( 72 Sr). 
Equation (3.18) with ( |3.19| ) and (3.20) proves the lemma. 



(3.20) 
D 



4 Operator-valued Riemann— Hilbert problem 

In this section we establish the operator-valued Riemann-Hilbert problem associated with the 
generating functional ( |1.3[ ). Let us consider the 4x4 matrix integral operator x ( m '(z\s,t) that 
has the following properties. Hereafter not only the index m but also variables s, t are omitted 
if unnecessary. 

1. x( z ) is analytic for z G C\C. 

2. Let Xext( z ) an d XintC- 2 ) De the limits of x( z ) horn outside and inside of the unit circle, 
respectively. It then follows that 

X ex t(z)=Xint(z)L(z) (zgC) (4.1) 

where the conjugation matrix L(z) is given by 

L(z) = l-l(z) l jk (z\s,t) = aj (z\s)b k (z\t). (j,fc = l,...,4) (4.2) 

3. x(oo) = 1. 

We point out that the conjugation matrix L(z) is the 4x4 matrix integral operator with variables 
s,t. For example, the condition (4.1) means 

(Xcxt) jk(z\s,t) = / drY](Xint)jn(z\s,r)L nk (z\r,t). (j, k = 1,... ,4) (4.3) 

Jo ~, 



n=l 



The connection of the operator- valued Riemann-Hilbert problem 1-3 to Lemmas ||-g in the 
previous section is summarized in the following theorem: 

Theorem. Suppose that the solution of the Riemann-Hilbert problem 1-3 exists and is unique. 
Define new vectors by 

roo /'OO 

a(z\s)= dtx(z\s,t)a(z\t) f3(z\s) = dtb(z\t)x~ l {z\t,s). (4.4) 

Jo Jo 

and introduce the integral operators A(s,t) and B(s,t) as Definition^ again. Then they satisfy 
Lemmas y 0. 

Proof. We start from the proof of Lemma pi. Using the definition of A(s,t) and the canonical 
integral representation of x{ z ) ( see IE % H) 

we obtain 



Since I (z) = (which is clear from the integrable condition (|2.7|)) there exists the inverse of 
the conjugation matrix: L~ (z) = 1 + l(z). The canonical integral representation of X~ l { z ) ls 
thus given by 

*" 1 <°)- 1+ sf/ t ,w «' i w- 1+ sf b - (4 ' 7> 



The identity x(0)% _1 (0) = X _1 (0)x(0) = 1 corresponds to Lemma ||. 

Let us compute the integro-difference relation of X ( z ) with respect to m and prove 
Lemma ^. We introduce the following integral operator 

^ m \z)= X {m) {z){ 11 +z- m l2 ). (4.8) 

This obeys 

*S(z) = *l™\z)L (z) (zeC) (4.9) 

with the conjugation matrix 

L (z) = (71 + z m l2 )L^(z)( 11 + z- m l2 ). (4.10) 

From the definitions of a^ m '(z), b^ m '(z) and L^ m '(z), we see that Lq(z) is independent of m. 
Hence applying Liouville's theorem we have 

^(™+i)( z )(^M( z ))-i = ll + i x (m+1) (0)7 2 (x (m) (0))- 1 (4.11) 

z 

which implies 

X {m+1) (z)(^i + -^72) = (v^7i + ^X (m+1) (0)72(x (m) (0))- 1 )% (m) (^). (4.12) 

\/x Jx 



By letting aS m >(z\s) act on this relation from the right the integro-difference equation ( |3.lD is 
derived. Similarly ( |3.2|) is obtained. Thus Lemma || has been proved. 

We show Lemma ||. In the same way as the FSFP case || and the impenetrable Bose gas 
case [10], in the neighborhood of C, the integral operator x( z ) can be decomposed into 

X(z) = x(z)Xo(z) (4.13) 

where x(z) is a single-valued, invertible and analytic in the neighborhood of C. Xo( z ) ^ s repre- 
sented by 

Xo(z) = l + ^log^Pl(z). (4.14) 

Since I (z) = its logarithm-derivative is computed as 

Due to this relation and Liouville's theorem the logarithm-derivative of x( z ) is written as 

d^x{z)x~\z) = £ ~^*e (4.16) 



£=± 

with the coefficient X±: 



Z-Ce 



X± = lim (z-C±)x(z)d i) Xo(z)Xo 1 (z)x \z) 

= -^x((±)i((±)x-\{±) 

= -^(CiWCi)*- 1 ^)- (4.17) 

10 



By definitions of a(z\s) and j3{z\s) the elements of X± are expressed by 



< 

2tt 



(X ± ) jk (s,t) = -^a j ((±\s)(3 k (( ± \t). (j,fc = l,...,4) (4.18) 



Let x( z ) ac t on ( 4.16| ) from the right and use this representation of X±. The differential 



equation ( |3,9| ) is thus obtained. Similarly ( p. 10 ) can be derived. The proof of the theorem has 



been completed. □ 

In consequence of Riemann-Hilbert problem 1-3 Lemmas |lj, [5| certainly hold. On the basis 
of these lemmas we can evaluate the large m asymptotic behaviors of correlation functions for 
the XX Z spin chain. 

5 Concluding remarks 

In the present paper we have derived a system of integro-difference equations that prescribe the 
generating functional of correlation functions and have established an operator- valued Riemann- 
Hilbert problem for the XXZ spin chain. It can be easily checked that, in the limit a — > — oo, 
our problem reduces to the problem associated with the FSFP obtained in ||. We are in a 
position to evaluate the long-distance asymptotic behavior of any correlation function for the 
XXZ spin chain. In a forthcoming publication our Riemann-Hilbert problem will be used to 
compute the large m asymptotic form of the generating functional. 

In the free-fermionic point r\ = 37r/4 the XXZ spin chain is reduced to the XXO model. The 
FSFP of the XXO model is known to be connected with an Riemann-Hilbert problem whose 



conjugation matrix is not an integral operator [11]. Its large m asymptotic form is computed as 

follows: 

/h + 2\ m2/2 
Pxxo(m) ~ [-^) (5-1) 

where h is the magnetic field. It is interesting to reproduce the same result from our Riemann- 
Hilbert problem by taking the limit a — ► — oo and r\ — ► 37r/4. Furthermore, we can evaluate any 
correlation function for the XXO model. 

Operator-valued Riemann-Hilbert problems appear in several subjects, apart from the com- 
putation of correlation functions of quantum solvable models. For example, it is well known that 
the classical inverse problem for integrable equations in 2 + 1-dimensions (Davey-Stewartson, 
KP, etc) can be expressed as an operator-valued Riemann-Hilbert problem. The investigation of 
operator-valued Riemann-Hilbert problems is thus important and interesting from standpoints 
of not only physics but also mathematics. 
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